Particle creation in a Robertson-Walker Universe revisited 
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We reanalyze the problem of particle creation in a 3+1 spatially closed Robertson- Walker space-time. We 
compute the total number of particles produced by this non-stationary gravitational background as well as the 
corresponding total energy and find a slight discrepancy between our results and those recently obtained in the 
literature Ql. 
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I. INTRODUCTION 

In a curved space-time the concept of particles is a subtle 
one and the meaning of what is a particle or what is a detec- 
tor becomes much more difficult than in a flat space-time J2l 
. Basically, this happens because in a general curved space- 
time we do not have Lorentz symmetry anymore and it is pre- 
cisely this symmetry that, in a flat space-time, allows us to 
identify the best vacuum state for the theory. Lorentz sym- 
metry assures that the vacuum state is the same for all inertial 
observers. 

Fortunately, in many problems of our interest, there are re- 
gions where a choice of the physical vacuum is quite natural. 
That happens when the time dependence of the field in these 
regions is harmonic (or, at least, almost harmonic). A basic 
condition to have a field with a harmonic time dependence is 
that the background metric and the eventual boundaries do not 
depend on time. In some specific cases, the background met- 
ric respects this condition only in asymptotic times (distant 
future and remote past). In these cases, the physical inter- 
pretation of particles is possible at those asymptotical times, 
but not between them. As a consequence, a non-static curved 
space-time background may lead to the phenomenon of parti- 
cle creation. 

The first one to discuss the problem of particle creation 
due to a curved space-time background was Schrodinger in 
1939JH, but the first one that carefully investigated this phe- 
nomenon was Parker in the late 60s 10]. The particle cre- 
ation in a 1 + 1 spatially closed Robertson- Walker space-time 



was investigated in J3, ID]. The 3 + 1 version of this prob- 
lem has been investigated in a recent work |Q]] @]. However, 
when generalizing the 1+1 solution to the 3+1 one, this author 
missed a degeneracy factor leading to an incomplete answer 
for the total number of particles created as well as the corre- 
sponding total energy. 

The main purpose of this short paper is to compute the to- 
tal number of particles created and the total energy related to 
them taking into account the correct degeneracy factor that 
arises when we are in a 3+1 dimensions. With the aid of ap- 
propriated graphics, we also make a brief discussion on how 
the parameters that appear in the metric affect the total number 
of particles produced. 



H. PARTICLE CREATION 

We consider the case of an expanding spatially closed 
Robertson- Walker universe whose line element and scalar cur- 
vature are given, respectively, by 

ds 2 = a 2 {^(drf - dp 2 - sen 2 p{d9 2 + sen 2 9 dtp)}, 
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where a(rj) is the scale factor, 77 is the conformal time and 
< p < 7r. In this case the equation of a massive scalar field 
conformally coupled to the metric is written as 



<9 2 itk 2 da duk (22 1 ^ 2fl 



t ~JT^ h m a + 1 + u k - V ang u k = 0, (2) 

arf a orj orj \ a orj z J 

where V arig is the angular part of the Laplacian operator on a 3-sphere. This operator has the hiper-spherical harmonics as 
eigenfunctions, that satisfy the following differential equation J3l 

Va ng Y(7, m 1 ,m 2 ; <p, 9, p) = -1(1 + 2)Y(Z, mi,m 2 ; <p, 9, p), (3) 

I 

where I = 0, 1, 2, mi = 0, ±1, ±2, ±Z and m 2 = gests the ansatz, 

0, ±1, ±2, ±mi. Hence, the V ans in equation (|2| sug- ]_ 

u k {ip, 0, p, 77) = — r^Y(l, mi,m 2 ; tp, 9, p)gi(rf). (4) 
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Substituting in (O, we obtain the differential equation for 



d 2 m(v) 

drf- 



where 



(5) 



(6) 



Note that this equation is similar to that of a mechanical har- 
monic oscillator with a time-dependent frequency. 

Now, we consider an exactly solvable case due to a conve- 
nient choice of a(i]), namely 



a(rj) 



'A + Btanh 



(7) 



where A, B and 7/0 are constants and A > B. Since we are 
considering an expanding universe, B > and r)o > 0. An 
inspection in FigureQ]allow us to interpret 770 as the time scale 
of the expansion of the Universe while A — B and A + B are, 
respectively, the scales of the size of the universe before and 
after the expansion. 




FIG. 1: Graphic of a 2 (77) versus r\ for a closed expanding universe. 



As we can see from Figure Q] there are two asymptotic 
times where the background metric is almost static, namely, 
the remote pass, characterized by 77 <C —770 and the distant 
future, characterized by 77 770. At those asymptotic times 
the corresponding vacuum states may be well defined if the 
solutions of <(5j have the following asymptotic limit 



Urn gj*<%) 



u! 



(PJ) 
I 



= v/(Z + l) 2 + m 2 (ATB), (8) 



where the superscripts p and / stand for past and future solutions, respectively. The solutions of (O that respect the asymptotic 
limits © are given Respectively, by 



(2^fiV k ) 1 



AT = — 1/2 2Fl I 1 + lLU i Vo, ^ W, 1 - itfVo; 1 - , 



ffi = ; -775 2-Fi (1 + i^i Vo,i^i Vo; 1 + iu{ i?o; £) , (9) 



where 2-F1 is the hipergeometric function |Q] and we have defined uj[^ := \ (oj( ± u}f\ and £ := (l — e 2, 'A'") 1 . Since u£ 

and u£ are not equal, the corresponding Bogolubov coefficients must be non-vanishing. Using some well known proprieties of 
the hypergeometric function, we can write u£ in terms of u£ and u£* as follows 

= y~] (akk'W k / + /3kk'W k *j , (10) 



k' 



with the Bogolugov coefficients given by 



Qk'k 



I uj( r(l - imf7j )T(-iuj{rio) ^ 



R r \<4 T(1-uJ?tio)T(uj{tio) 

Pk'k - OkkM/ —p ,, _ — . , . _ — r- (ii) 
V w ; r(iw z 77o)r(l + iu; ; 770) 
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In the remote past all the inertial particle detectors register the complete absence of particles in the state |(P) (the vacuum state 
associated to However, in the distant future, any inertial particle detector will register a number of particles with quantum 
numbers k in the |0 P ) state given by 



(w\nI\w>) = J2\M 2 = 



sinh 2 (7rtj[ '770) 
sinh(7ra;| > 77o) sinh(7rw^?7o) 

The total number of produced particles and the total energy associated to them are given, respectively, by 

sinh 2 (7rwj ( ~ ) » 



(12) 



(0*|iV/|OP>=££|/3 k v| 2 = £0 + l) : 



z=o 



Vo) 



sinh(7ro;['77o) smh^Trujf rjo) 



(13) 



e = Yw\nlW)<4 



k 



r 



, Lof sinh {ttojI 770) 
sinh(7rwf 770) sinh(7ro;/r7o) 



(14) 



Since we were not able to perform summation dT3b analyti- 
cally, let us make a numerical analysis of the preceding results 
constructing, for instance, the graphic of the total number of 
particles created TV versus the mass of the field excitations m. 
Naively, we could expect that the total number of particles cre- 
ated was a monotonically decreasing function of m. However, 
as we can see from the Figure |2j there is a value mo for the 
mass of the field at which the total number of created particles 
N(nio) reaches a maximum. 




mo 



FIG. 2: Graphic of the total number of particles produced versus the 
mass related to the field. 

It is convenient to define a width Ajtio as it is done in a 
With this purpose, we define m + and 



Gaussian distribution. 
m_ such that 



N(mo + m+) = N(mo — m_) 



l N(m ), 



(15) 



so that Amo := ?n+ — m_. 

Note from equationdDb that N(mo), mo and Amo depend 
only on A — B, A + B and 770- Despite the fact that we were 
not able to get a closed analytic expression for mo, Amo an d 
TV (mo) as a function of the metric parameters A — B, A + B 
and 7]q, we can see numerically how these quantities depend 



on those parameters. We do that by plotting the graph of the 
total number of created particles as a function of the field mass 
for several values of A — B, A + B and 770 (Figure[3]). At first 
we plot a control curve (in Figure [3] it is represented as a 
continuous line), and then we modify one of the parameters, 
plot a new curve and compare to the control one. We repeat 
this procedure for all parameters. 
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FIG 3: Graphic of the total number of created particles N versus 
the field mass m for different values of parameters A — B, A + B 
and r)o . The continuous line represents the case where rjo = 0.01, 
A — B = 1 and A + B = 101. The dashed line represents the 
case where r; = 0.01, A - B = 3.25 and A + B = 101. The 
dotted line represents the case where 7/0 = 0.01, A — B = 1 and 
A + B = 51. Finally, the continuous dotted line represents the case 
where r/o = 0.0115, A-B = l and A + B = 101. 



Changing the parameters A — B, A + B and 770 as we did 
in Figure [3] and then analyzing numerically what happens to 
the graphic of N versus m, suggests the following behaviors: 
iV(mo) gets larger and A7« g ets smaller as A — B and r] 
decrease or A + B increases. On the other hand m gets larger 
as A — B, A + B or 770 increase. 
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m. CONCLUSIONS 

In this work we computed the total number of particles cre- 
ated and the total energy related to them in a 3 + 1 spatially 
closed Robertson- Walker space-time. By using appropriated 
graphics, we also discussed how the parameters that appear 
in the metric affect the total number of particles produced. 
We found an unexpected behavior of N as a function of m, 
namely, starting from zero, it increases until it reaches a max- 
imum value at mo and after that it decreases monotonically as 
m increases. 

The essential difference between our calculations and that 
presented in 0]] is that this author does not take into account 



the degeneracy factor (I + l) 2 that appears in equations < TT~3b 
and (fT4l . This factor must be included due to the degeneracy 
in the quantum numbers mi and m<x- We think our discussion 
may be of some help in the analysis of similar problems, as for 
example, if one considers other functions a(r\) with behaviors 
slightly different from that one considered here. 
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